Unbiased acceleration measurements with an electrostatic 
accelerometer on a rotating platform 

Benjamin Lenoir a , Bruno Christophe a , Serge Reynaud b 

a Onera - The French Aerospace Lab, 29 avenue de la Division Leclerc, F-92322 Chatillon, France 
b Laboratoire Kastler Brossel (LKB), ENS, UPMC, CNRS, Campus Jussieu, F-75252 Paris Cedex 05, France 



o 

(N 
C 



O 



(N 
> 
On 

On 
O 



X 



Published in Advances in Space Research 51 (2012) 188-197 
doi: |10.1016/j.asr.2012.08.012l 

6 January 2013 



Abstract 

The Gravity Advanced Package is an instrument composed of an electrostatic accelerometer called 
MicroSTAR and a rotating platform called Bias Rejection System. It aims at measuring with no bias 
the non-gravitational acceleration of a spacecraft. It is envisioned to be embarked on an interplanetary 
spacecraft as a tool to test the laws of gravitation. 

MicroSTAR is based on Onera's experience and inherits in orbit technology. The addition of 
the rotating platform is a technological upgrade which allows using an electrostatic accelerometer to 
make measurements at low frequencies with no bias. To do so, the Bias Rejection System rotates 
MicroSTAR such that the signal of interest is separated from the bias of the instrument in the 
frequency domain. Making these unbiased low-frequency measurements requires post-processing the 
data. The signal processing technique developed for this purpose is the focus of this article. It allows 
giving the conditions under which the bias is completely removed from the signal of interest. And 
the precision of the unbiased measurements can be fully characterized: given the characteristics of 
the subsystems, it is possible to reach a precision of 1 pm s -2 on the non-gravitational acceleration 
for an integration time of 3 h. 

Keywords Electrostatic accelerometer; Rotating platform; Bias rejection; Absolute measurement; 
Modulation 

PACS 02.50.Ey; 04.80.Cc; 06.30.Gv; 07.87,+v 

1 Introduction 

The experimental tests of gravita tion are in good agreement with its current theoretical formulation 
referred to as General Relativity (jWi But contrary to the quantum description of the three 

other fundamental interactions, it is a classical theory, which suggests that another description of grav- 
itation lies beyond General Relativity. From the experimen tal point of vieWj there are still open win- 
dows for deviations f r om G eneral Relativity at short range (jAdelberger et al. , 2003) and at long range 
( Jaekel and Revnaudl . 2005). Galactic and cosmic observations also challenge General Relativity. The 
rotation curves of galaxies and the relation between redshifts and luminosities of supernovae, which 
are interpreted as m anifestations of "dark matter" and "dark energy" respectively ( Frieman et all [2008 : 
Bertone et al. , 20051 ) , may also be s een as a hint that General Relativity could be an imperfect description 
of description at these large scales (jAguirre et all l2001t iNoiiri and OdintsovL[2007h . 

In this context, testing General Relativity at the largest possible scales is essential. For man-made 
instruments, the Solar System can be used as a laboratory for gravitational experiments. NASA performed 
such a test of gravitat ion with the Pioneer 10 and 11 missions. The outcom e was a signal now known as 
the Pioneer anomaly ( Anderson et al. . 19981 l2002bt iTurvshev et al. . 2011). The long term variations of 
this signal may be explained as an anisotropic thermal effect ( Bertolami et al. , 2010t Rievers et al. , 20101 : 
Rievers and Lammerzahl . 2011 : Turvshev et al. . 2012) but periodic anomalies have also been identified 
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( Lew et al. . 2009t Courtv et al. , 2010f). Nevertheless, the Roadmap for Funda mental Physics in Space 
issued by ESA in 2010 ( Fundamental Physics Roadmap Advisory Team , 2010| ) stresses the importance 
of testing gravitation at large scales with missions to the outer planets. To do so, it recommends the 

development of accelerometers compatible with spacecraft tracking at the 10 pm s~ 2 level . 

Se veral missions have already been proposed ( Anderson et all 2002at Dittus et al. , 20051 : Johann et al 



20081: iBertolami and Paramosl 120071 : IChristophe et all l2009t IWolf et all 120091) with the aim of improv 
ing the knowledge of the gravitational field in the Solar System. Many of them propose to embark 
an accelerometer which will measure the non-gravitational forces acting on the spacecraft in order to 
distinguish unambiguously the non-gravitational accelerations from gravitational effects. In practice, 
acceler ometers measu r e a combin ation of the spacecraft non-gravitational acceleration and additional 
terms (C arbone et al. 1. 120051 120071) . It will be assumed in this article that these additional terms either 
are negligible or can be corre cted, such that the external signal will be referred to as the spacecraft 
non-gravitational acceleration ( Lenoir et al. . 2011bl ). 



This ar t icle d eals with the Gravity Advanced Package, an instrument proposed on Laplace mission 
(|BiesbroekL l2008h . which is designed to make measurements of the non-gravitational acceleration of a 
spacecraft with no bias. It provides an additional observable which measures the departure of the space- 
craft from geodesic motion. To do so, post-processing is required, which is the subject of this article. It 
allows retrieving separately the acceleration without bias and the bias of the instrument, these quantities 
being referred to as "post-processed" quantities. 

In a first part, the Gravity Advanced Package will be presented as well as its performances and 
the measurement principle. Then the post-processing method will be described and conditions will be 
given so that the bias can be effectively removed from the measurement. These conditions will allow 
deriving measurement procedures to make unbiased measurements. The emphasis will then be put on the 
characterization of the post-processed quantities, i.e. the quantities without bias, and on their uncertainty. 

This paper is focused on the performances of the accelerometer. The constraints in the integration of 
the instrument in the spacecraft with the a im of preserving these performances requires additional work. 



The OSS mission (jChristophe et all 120121 ) proposes a spacecraft design which takes into account these 



concerns. 



2 Instrument principle, design and performance 

The Gravity Advanced Package is made of an electrostatic accelerometer, called MicroSTAR, which can be 
rotated with the Bias Rejection System. This technological upgrade allows removing the bias introduced 
by MicroSTAR. 



2.1 Overview of the instrument 

MicroST AR is a 3-axes electrostatic accel erom eter (Josselin et al. , 1999) based on O NERA's expertise in 
this field (|Touboul et all Il999t iHudson et all , 12007k iTouboul and RodriguesL l200lh . In orbit technology 
(CHAMP, GRACE and GOCE missions) is used with improvements to reduce power consumption, size 
and mass. 

The core of the accelerometer is composed of a proof mass inside a cage made of six identical plates. 
The motion of the proof mass with respect to the cage is detected by capacitive measurement. A control 
loop adjusts the potentials of the electrodes in order to keep the proof mass at the center of the cage. 
The numerical values of these potentials, which are the outputs of the instrument, are proportional 
to the components of the acceleration of the proof mass with respect to the cage on each axis of the 
accelerometer. 

The Bias Rejection System is a rotating platform composed of a rotating actuator and a high resolution 
angular encoder working in closed loop operation. Piezo-electric technology is envisioned for the actuator. 
It allows designing a device with no need for any kind of gear, so as to reduce mass and volume. The 
piezo-electric motor is operated in a slip-stick mode. Finally, even if piezo-electric motors have a non-zero 
torque in the power-off mode, a blocking system will be implemented to prevent unwanted motion during 
launch and maneuvers. 
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2.2 Performance of the electrostatic accelerometer 

The perform ance of MicroSTAR is measured via the power spectrum density of the noise on the measured 
acceleration ( Lenoir et al. . 2011hl Fig. 4). The analytic formula S n (in m 2 s~ 4 Hz -1 ) as a function of 



frequency, for a measurement range equal to 1.8 x 10 4 ms 2 , is 




/ \ , ( f 



with K = 5.7 x 10" 11 m.s- 2 .Hz~ 1/2 . 

In addition, the instrument has a bias. It corresponds to the deterministic low-frequency variations 
of the accuracy of the accelerometer. It is due to the gold wires which are used to keep the polarization 
of the proof mass constant and to the geometrical imperfections of the instrument. In previous missions 
relying on electrostatic accelerometers (CHAMP, GRACE and GOCE), this bias was not a problem since 
the measurement bandwidth was 0.1-100 mHz. On the contrary, for the application foreseen in this 
article, very low-frequencies measurements are to be made and it is therefore required to remove the bias 
from the measurements. 



2.3 Measurement principle 



The measurements made by MicroSTAR along its three axes x, y and z, which are supposed to be 
orthogonal, are 
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where Ski K , k 2K , b K and n K (k £ {x;y;z}) are respectively the scale factors, the quadratic factors, the 
bias and the noise on each axis. The bias has a deterministic temporal variation whereas the noise is 
a null-mean stationary stochastic process whose PSD is given by Eq. (fTJ). The quantities a K are the 
components of the non-gravitational acceleration in the reference frame of MicroSTAR. As far as orbit 
reconstruction is concerned, these quantities are not the ones of interest since the Bias Rejection System 
rotates MicroSTAR with respect to the spacecraft. The spacecraft is supposed to be stabilized along the 
three axes. The transformation matrix P moves a vector from the spacecraft reference frame (whose axes 
are X, Y and Z) to the accelerometer reference frame. In its simplest form (but without any loss of 
generality), the expression of P is 

[10 

P= cos(0) sin(0) (3) 
- sin(0) cos(0) 

where is a monitored angle, which measures the rotation of the accelerometer with respect to the 
spacecraft. Considering only the plane perpendicular to the axis of rotation of the accelerometer, Eq. (ffi) 
becomes 



m y = (1 + Ski y ) [cos(0)ay + sin(0)a^] 
+k 2y [cos(0)ay + sin(0)az] 2 + b y + n y 
m z = (1 + Ski z ) [— sin(0)ay + cos(0)a.z] 
+k 2z [- sm(8)a Y + cos(8)a z } 2 + b z + n z 



(4a) 
(4b) 



The measurements on the axes y and z are combinations of the quantities ay and az- These are 
the quantities needed so as to measure the impact of non-gravitational forces on the trajectory of the 
spacecraft. This fact associated with the possibility to give the angle any possible time variation allows 
measuring ay and az without bias. On the contrary, on the axis x, there is no possibility with this 
instrument to remove the bias from the measurement m x so as to retrieve ax- To do so, another rotating 
platform would be required. It is not the topic of this article but the method developed here can be 
applied to this more complex setup. In practice, the axes Y and Z will be in the orbit plane, in which 
non-gravitational forces are expected to impact the trajectory of the spacecraft. 



lr The orthogonality of the measurement axes depends on the orthogonality of the proof mass faces. Assuming a non- 
gravitational acceleration equal to 10~ 7 m s~ 2 (cf. section|3]2jl, the orthogonality needs to be controlled at a level of 10 uxad 
in order to have an accuracy on the measurement of 1 pm s — 2 . This level is technologically reachable. 
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In the following, it will be assumed that N measurements are made with a sampling frequency called 
f s . It corresponds to a time step called St — l// s . The scale and quadratic factors will be supposed to 
be constant. In Eq. (jj]), there are therefore AN unknowns (ay, az, b y , b z at each sampling time) and 
2N measurements {m y ,m z ) spoiled by noise (n y ,n z ). In the rest of this article, for each of these eight 
quantities as well as for 9, the notation x will be a vector of jMjv,i(K) whose components are the values 
of x at each sampling time and x k is the value of x at the sampling time k x St. 



3 Signal processing method 

The relations between the measurements made by the instrument, m y and m z , and the non-gravitational 
acceleration in the spacecraft reference frame, ay and az, has been expressed. In this section, the data 
processing method is presented. In particular, conditions are derived in order to remove the bias from 
the measurements. 

3.1 Linearization of the problem 

A first step is to linearize Eq. such that they can be written in a matrix form. To do so, it is assumed 
that k2 y — k2 Z = 0. This hypothesis will be shown in paragraph s. 1 1 not to be restrictive in the framework 
presented here. The two following diagonal matrices, belonging to Af^M), 



A c = diag[cos(6» fe )] and A s 
allow writing Eq. @ in the matrix form 



diag[sin(0 fc )], k& \\1;N\\ 
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(5) 



(6) 



(7) 



(8) 



The set of solutions for this system is infinite. It is the affine space X p + ker(J), where X p is a given 
solution of the linear equation. This formal resolution gives no useful information on the non-gravitational 
acceleration since it provides an infinite number of solutions. 



3.2 Generalized noise 

Before going further, it is necessary to consider the matrix J more carefully. When solving Eq. ©, J is 
supposed to be perfectly known. It is however not the case since the knowledge of the angle involved 
in the definition of J may suffer a bias and noise. There is a discrepancy between the true value of the 
rotation angle 9* and the measured one 9: 



59 



(9) 



where be is a bias and 59 a random process (whose mean value is equal to zero). Using the same notations, 
this leads to a noise described by 5A C and 5A S on the matrices A c and aJ3 

The impact of the bias on the precision of the measurement has been assessed in ( Lenoir et al. , 201 ltJ ) 
and it has been shown that be must be smaller than 10~ 5 rad in order to meet the expected performances. 

In order to take into account the impact of the noise 59. it is possible to introduce a generalized noise: 
the quantities n y and n z in equation © are replaced by n y = n y + n y and n z = n z + n z with 



n y = (1 + 5ki y ) [<5A c a Y + 5A s a z ] 
n z = (1 + Sku) [-<SA s a Y + 5A c sl Z ] 

2 SA C = diag[cos(£>k + 89^) — cos(0j.)] diag[— 89^ sin(6fe)] and <5A S = diag[sin(#k - 



(10a) 
(10b) 

e k )-sm(9 k )] «diag[50 fc co S (0 fe )] 
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This additional noise depends on the non-gravitational accelerations ay and az and on 5A C and 5A S . 
As a result, the smaller the magnitude of the external acceleration is, the smaller the noise due to the 
uncertainty on 9 is. 

It can be used to derive the requirements on 59 such that the predominant source of uncertainty is 
MicroSTAR and not the Bias Rejections System. To have such a result, one needs Sf l (f) -C S n (f) around 
the modulation frequency 1/t (cf. section f4 . 1 [) . where Sh(f) is the power spectrum density of the noise 
n due to the rotating platform (cf. Eq. (jTUJ) ) . Assuming that a y w a z s=s ajvGi we have h 59clng- This 
leads to the following requirement : 



V/e 



2r' 2~r~ 



Sn(f) 



"NG 



(11) 



To compute a^G, it is assumed that the main contributor is solar radiation pressure and that the space- 
craft is at one astronomical unit (called do) from the Sun. The power car ried by solar photons by surfac e 



unit at this distance is approximately equal to P — 1.366 x 10 3 W m~ 2 (jWillson and Mordvinovl l2003f l. 
Conside ring a ballistic co efficient equals to Cb = 0.1 m 2 kg -1 , which is the order of magnitude for Laplace 
mission (jBiesbroek , 2008), the non-gravitational acceleration is equal to clnq — CbP/c = 4.6x10 7 m s 2 
at one astronomical unit, where c is the speed of light. Taking the minimum value of S n , the requirement 
on S$e reads: 

(12) 



y/Sse(f) < 1-3 x 10~ 4 rad.Hz 



-1/2 



In the rest of the article, it will be assumed that this condition is verified and only the noise of MicroSTAR 
will be considered. 



3.3 Conditions for bias rejection 

The general approach presented above to solve the linear system does not give useful information on the 
non-gravitational acceleration or on the bias of the instrument. Since it is impossible to obtain the value 
of the unknown quantities at each sampling time, it is necessary to narrow the information retrieved from 
the data. A possibility is to look for the projection of the vectors ay and az on a vector subspace (of 
dimension p a < N) whose basis is made of the column of a matrix V a € Ai M,p a (K) , which are supposed to 
be orthogonal for the usual scalar product on Mr. As a result, the goal is to find the numerical values of 
V' a ay and V^'az knowing m y and m z [M 1 is the matrix transpose of M). In this article, the choices of V a 
will allow retrieving the mean value of the acceleration without bias and the slope of the acceleration over 
one modulation period. But other choices of V a can be made to retrieve for example sinusoidal variations 
of the signal. 

Under the following four conditions on the bias, the angle 9 and the projection matrix V a 

V a 'A v h K = 0, with v € {c; s} and k £ {y; z}, (13) 

and assuming that 8ki y = 5ki z — Ski, the unbiased values of the external signal can be recovered: 

f K'(l + MiW = K'A c m y - y Q 'A s m z (14a) 
I K'(l + Sh)a z = K'A s m y + K'A c m z (14b) 

Calling Vfe the A:-th column of V a , the conditions (|13[) can be expressed in the frequency domain 

(T St {v fc cos(0)}, F st {b K }) = (15) 

where J-$ t is the discrete time Fourier transform and (•) is the usual scalar product. This equation means 
that the bias and the modulated signal must be orthogonal in the frequency domain. 

It is a priori not possible to know whether conditions (|13p are fulfilled since the temporal evolution 
of the bias of the instrument is not controlled. However, as already mentioned, the bias corresponds 
deterministic low frequency variations. It is therefore possible to assume that b y and b z belongs to a 
vector subspace defined by the columns of V& € (R) {p < N). Given this hypothesis, conditions 

(fl3|) come down to 



f V a 'A c V b = 

\ V a 'A s V b = 



(16a) 
(16b) 
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The results presented in this section can be found by solving equation ^ with a modified least square 
method. The matrix J, which is unknown (because of the scale factors), is replaced by the matrix 

J G M-N,2{p a +p b ) 

' A c V a A s V a V b 0" 
-A s V a A c V a V b 



J 



(17) 



with V b G M N:Pb (M 
and that b y and b 5 



) . And it is assumed that ay and az belong to the subspace generated by V a G A4 M,p a 
belong to the subspace generated by Vb- Section [531 will build on this approach. 



4 Unbiased measurements of non-gravitational acceleration 

Based on the conditions (|16l) required for a correct demodulation and given some assumptions on the 
matrices V a and Vb, it is possible to design a calibration signal, i.e. a pattern for the angle 0, which allows 
for completely removing the bias from the measurements. 



4.1 Choice of a calibration signal 

The calibration signal looked for will be periodic, with a period called r. First, some practical concerns 
restrict the possible pattern. Because it can be assumed that rotating the accelerometer will induce 
vibrations and therefore spoil the measurements, the angle 9 will have to be constant when the measure- 
ments are done. As result, calibration signals such that 9(t) = 2wft, where / is an angular frequency, are 
forbidden. Moreover, because the accelerometer may not be perfectly centered on the rotating plate, the 
rotation induces Coriolis and Centrifugal forces which spoil the signal. Finally, rotating constantly may 
lead to a quicker breakdown of the instrument. 

Another practical concern, which appears if no slip ring is used, is about the wires between the ac- 
celerometer and the spacecraft. Because of them, it is not possible to rotate the accelerometer indefinitely. 
Therefore, the angle 9 will have to stay in the interval [0; 2ir}. 

To go further, it is necessary to be more specific on the matrices V a and Vb. First, constant values of 
the non-gravitational acceleration during each modulation period will be looked for and the bias of the 
instrument will be supposed to be, for each period, an affine function of temperature. Therefore, 



V a 







and Vb 







(18) 



where 1 is a matrix of Ai 



9,11 



whose coefficients are 1, and T is a matrix of M. 



N,l\ 



made of the 



values of the temperature at each sampling time. The integer q is the number of sampling points in one 
period. It is assumed that r and f s are such that rf s is an integer and q — rf s . In this approach, the 
variation of temperature will be assumed to be driven by the heat generated by the rotating platform: at 
each rotation, heat is generated and induces a temporary increase of temperature. 

Figure Q] shows two examples of calibration signals which fulfill the conditions (fH)]) under the previous 
assumptions for the bias, the non-gravitational acceleration and the temperature. Let us consider the 
signal of Fig. 1(a) and go back to the assumptions made previously on the linear and quadratic factors^. 
In Eq. the quadratic terms are constant because these two equalities are always true: sin(6>) = 
and [cos(#)] 2 = 1. Therefore, the quadratic terms behave as a bias which will be separated from the 
non-gravitational acceleration. Concerning the assumption on the equality of the scale factors, it has to 
be noticed that A c 2 = Id^r and A s — 0. Therefore, the derivation leading to Eq. (TH)) still hold without 



the assumption on the scale factors. On the contrary, for the signal of Fig. 1(b) as well as for any signal 
for which 9 has values different from ° and 180 ° , these remarks on the scale and quadratic factors do 
not apply. As a conclusion, only signals for which measurements are made when 9 = 0° and 9 = 180 ° 
should be considered. 

The previous hypothesis made on V a and Vb allows to derive simple calibration signals. They are 
however restrictive because it is assumed that during a modulation period the signal and the bias are 



3 It was assumed that k% y = k<2 z = and Ski y = 5ki z 
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Figure 1: Example of two calibration signals 6{t) which fulfill the conditions given by equations (fT6j) for 
the matrices V a and Vb defined by equations (fT5)h The rotating duration corresponds to 33.3 % of the 
modulation period r = 1 arbitrary unit. Two periods are represented, separated by circles (o). 
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Figure 2: Example of a calibration signals 6(t) which fulfill the conditions given by equations (fTH|) for 
the matrices V a and Vb defined by equations (fT^|) . The rotating duration corresponds to 33.3 % of the 
modulation period r = f arbitrary unit. Two periods are represented, separated by circles (o). This 
signal is different from the one of Fig. 1 (a) because it does not display a periodicity of 0.5 arbitrary unit. 



constant (with a temperature dependence for the bias). To go further, it is possible to design a calibration 
signal assumi ng that the bias is for each period an affine function of time but does not depend on 



temperature (|Lenoir et all l2011al) . With this signal the mean and the slope of the non-gravitational 



acceleration on each modulation period will be recovered. In this case, 



V a = V b 



t„ 



1, o 







(19) 



where t q is a matrix of At g ,i(M) such that t q , = (k — q/2)5t. Figure [2] shows a calibration signal which 
fulfill conditions (fT6|) . Note that the remarks made on the scale and quadratics factors hold for this 
calibration signal. Contrary to the calibration signals of Fig. [U the pattern in this case depends on the 
masking time which is introduced in the following section. In the rest of this article this calibration signal 
will be used. 

In case the bias of MicroSTAR does not belong to the subspace generated by Vb, then the signal of 
interest is not perfectly recovered: it is spoiled by the quantities V a 'A„b K (y € {c; s}, k £ {y; z}). 
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4.2 Masking 

As mentioned in the previous paragraph, measurements made when the accelerometer is rotating are not 
considered for data reduction because they may be spoiled by unwanted signals. Therefore, during post- 
processing, the data acquired when the accelerometer is rotating must not be taken into account. This 
will be refered to as "masking". To introduce this masking feature in the signal processing, let consider 
the diagonal matrix M € A^at(K) defined by: Mf-k = 1 if 0fc = and 9^ = 0, and M^h — otherwis^]. 
Then in Eq. (TT3|). the matrix V a is replaced by V a = MV a . 

The duration of masking is a key parameter in the precision of the post-processed quantities: the 
longer it is, the more data points are lost and the uncertainty increases (cf. section l5^f|) . The total 
duration of masking during one period is called Tm- 



5 Demodulated quantities 

The demodulation signals introduced in the previous section allows to retrieve unbiased measurements 
of the non-gravitational acceleration of the spacecraft. The focus will be now to characterize these post- 
processed quantities in term of uncertainty. 

5.1 Autocorrelation of the non-gravitational acceleration mean 

The calibration signal of Fig. fallows to recover affine variations of the external signal on each modulation 
period. In term of spacecraft navigation, the goal of the instrument is to measure the impact of non- 
gravitational forces on the dynamics of the spacecraft. And the variation of momentum during one 
modulation period Ap^vG of the spacecraft due to the non- gravitational forces Fng 1S equal to the mean 
of the non gravitational forces times the modulation period: 

ta+r 



& PNG = / F NG (t)dt = t (F NG ) (20) 

J t T 

where to is an arbitrary time and (■) is the mean during a duration r. 

As a result, only the mean values of the external signal are of interest, and the subsequent analysis will 
be restricted to the matrix V a defined by equation (|18p . Under the assumption introduced previously, the 
demodulated acceleration are defined by Eq. (|14[) . In order to have normalized quantities, it is necessary, 
as in the least square method, to multiply this equation on the left by (V^'V^) -1 . Under the assumptions 
considered here, this matrix is diagonal with all the coefficients equal to q = | | 2 , where S A^at,i(M) 
is ith column of the matrix V a . q is independent of the index i. 

Let us call Cj = A c Vi £ A^at,i(K) the ith column of the matrix A C U Q , ay^ = (1 + Skxy^^ay / q the ith 
component of the column vector (1 + 5fcij / )(V r a 'U a ) _1 VJ J 'aY, and ozi = (1 + <5fcl 2 )v^az/(j , the ith component 
of the column vector (1 + (5fci 2 )(V r a 'V r Q )^ 1 V r a 'a Z - 

The quantities 5y, and azi are the means of the non-gravitational acceleration of the spacecraft for 
the modulation period i along the axes Y and Z respectively. Under the assumption made earlier, the 
accuracy of the measurements is perfect, i.e. their expected values is equal to the true values. Concerning 
the precision, assuming that n y and n z are independent and have the same power spectrum density, S n 
defined by equation (fTJ), the covariances between the post-processed quantities are 



Cov(a Y i,a Yj ) = Co\(azi,a Zj ) = J ^ S n (f) I I d f ( 21 ) 

and 

Cov(2? i ,o^ J ) = (22) 

The result given by equation (|2~T|) is true only if the signal has been filtered before digitization by a 
perfect low-pass filter with a cut-off frequency of f s /2 so as to avoid aliasing. 

According to Fig. [31 the integral of equation |2T|) select the noise power spectrum density at the 
frequency 1/r and approximately integrate it on an bandwidth 1/r for i = j. In order to minimize the 
absolute value of the covariance, it is therefore necessary to select the noise at the frequencies where it is 



4 9 and 9 are respectively the angular velocity and the angular acceleration of the rotating platform. 
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Figure 3: Re Tstic^Tsticj} for the calibration signal of Fig. [5] with a modulation period t of 600 s, a 
masking time Tm of 200 s and a sampling frequency f s of 1 Hz: (a) i — j, |(b)| j = i + 1. For figure (a) 



the peak is approximately at the frequency 2/t and its frequency width is approximately 1/r. 



minimum, i.e. for / G [10~ 2 ; 2 x 10 _1 ] Hz, which correspond approximately to modulation period between 
5 s and 100 s. Too short modulation periods are impossible to implement in practice. Therefore, in the 
following, a modulation period equal to 10 min will be considered. 
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Figure 4: Autocorrelation function Rd[k] defined by equation (|23[) for the calibration signal of Fig. [5] with 
a modulation period r of 600 s, a masking time Tm of 200 s and a sampling frequency f s of 1 Hz. 



The main interest of the demodulation process is to know the mean acceleration over a modulation 
period r. This process gives birth to two new discrete-time quantities ayi and azi indexed formally by 
i g Z. It is possible to introduce the autocorrelation function Rd [k] which is the same for both quantities 
and which is defined by 

R d [k] = Cov(a^ i+k , avi) = Cov(a^ i+k , a^i). (23) 

Fig. [3] shows that the autocorrelation function is close to the one of a white noise. This means that 
the post-processed quantities are approximately independent. In term of power spectrum density, this 
corresponds to a level of 10 -10 m s~ 2 Hz -1 / 2 with a cut-off frequency equal to 8.3 x 10~ 4 Hz. Since the 
uncertainty on the demodulated accelerations is known and characterized, it is now possible to use them 
to gain more information on the non-gravitational accelerations. 
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Figure 5: Uncertainty on the mean acceleration as a function of the integration time assuming that the 
post-processed quantities are independent or have the autocorrelation function plotted in Fig. 01 The 
plot is for the calibration signal of Fig. [5] with a modulation period t of 600 s, a masking time Tm of 
200 s and a sampling frequency f s of 1 Hz. Its shows that the independence approximation is correct. 
For an integration time of 3 hours, the uncertainty is equal to 1 pm.s~ 2 on the mean acceleration. 



5.2 Further characterization of the non-gravitational acceleration 

In order to increase the precision, it may be interesting to know the mean acceleration over periods of 
time longer than the modulation period. To do so, one needs to average the demodulated accelerations 
over the period of time of interest. Figure [5] shows the uncertainty on the mean acceleration for different 
integration time. As the noise is nearly white, the uncertainty on the mean decreases as 1/yT where T 
is the integration time. 

It is also possible to look for sinusoidal variations of the non gravitational acceleration with a known 
frequency /*. The goal is to find the coefficients a and j3 of the time varying signal a cos(2ir f*t) + 
f3 sin(27r/*t) using the values of the non-gravitational acceleration for each modulation period. According 
to the Nyquist-Shannon theorem, it is not possible to recover sinusoidal variations at frequencies higher 
than half the sampling frequency, i.e. /* > l/(2r). Conversely, when l/(r/*) becomes too large, the 
uncertainty diverges. The value for which this happens depends on the number of post-processed points 
used to fit the sinusoidal variation: the more points are used, the easier it is to fit low frequency signals. 
For r = 10 min and /* the frequency related to the revolution period of the EartrH, the frequency ratio is 
l/(r/*) = 144. In this particular configuration, one obtains with 60 points and for a modulation period 
of 10 minutes (which corresponds to 10 hours of measurement): 

!VCov(a, a) = 8.7 x 10" 13 m.s~ 2 (24a) 
y/Cov{P,P) = 7.3 x 10" 13 m.s~ 2 (24b) 
v/-Cov(a,/3) = 2.4 x 10~ 13 m.s~ 2 (24c) 

These values show that it is possible to obtain, in this configuration, the amplitude of the sinusoid with 
a precision better than 1 pm s~ 2 . 

5.3 Generalized least square/optimal filtering 

In section 13.31 it was mentioned that the process described until now corresponds to a least square 
(LS) method. This method provides estimates with a minimum variance only when the noise is white. 
However, the noise of MicroSTAR does not fall in this category. The generalized least square (GLS) 

5 Because the instrument will be used for spacecraft tracking, one may be interesting in sinusoidal variations at the 
revolution period of the Earth. These variations are also of interest for the fundamental physics objectives discussed in the 
introduction. 
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method (|Cornillon and Matzner-Loberl 120071 ) provides an estimate with minimum variance whatever the 
measurement noise is. 

This method is similar to the optimal filtering technique (jPapouhsl . 119771 p. 325): the first one is 
expressed in the time domain whereas the second one is express in the frequency domain. It is possible 
to express the components of the inverse of the covariance matrix Vqls = (J'^ 1 J) 1 using the power 
spectrum density of the noise instead of its covariance matrix: if v and w are two column vectors, then 



w = 



-—T St { v }(f)F 5t {v,}(f)df 

-3k 



(25) 



One may process the data form the accelerometer using the generalized least square method. However, 
in the specific case of the problem studied here, the gain is rather small. Indeed, Fig. [3] shows that, for 
the calibration signal considered, the Discrete Time Fourier Transform is peaked around the frequency 
2/r and the noise PSD does not vary much on the interval [10~ 3 ; 10 _1 ] Hz. As a result, using Parseval 
theorem, 




' \F St { Ci }{f)\ 2 df 



1 



Sn(f) 



Sn(f) 



•MciX/) 



q.St 



df 



(26) 



Therefore, the autocorrelation function plotted in Fig. 2] is nearly the same as the one obtained with 
the GLS approach: the autocorrelation function obtained with the GLS approach is similar to the one 



of a Gaussian noise and its value for k = is 1.71 x 10 



-23 



s" 4 instead of 1.84 x 10" 



for 



Fig. 01 The difference in the level of precision on the post-processed quantities is also visible on Fig. El in 
section 15.41 



5.4 Optimization of the masking time and calibration period 





Modulation period (min) 
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Modulation period (min) 
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Figure 6: Uncertainty on the demodulated acceleration for an integration time of one hour and for the 
calibration signal of Fig. [5] as a function of the modulation period. The plots are parametrized by the 
masking time. The dash line shows the modulation period which gives the minimum uncertainty for 
each masking time. Figure (a) is obtained using the least squares method to process the data whereas 



figure [(b)] is obtained with the generalized least squares method. 



Until now, only one modulation period (r = 10 min) and one masking time (Tm = 200 s) have been 
considered. But since their value impact the uncertainty on the demodulated accelerations, it is legitimate 
to choose these values such that the uncertainty is minimized. 

Figure [H] shows the uncertainty on the demodulated acceleration for an integration time of one hour 
and for the LS (a) and GLS (b) methods. This value is computed by taking the numerical value of Rd[0] 
and by multiplying it by y^r/T, where r is the modulation period and T = 1 hour. This assumes that the 
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demodulated accelerations are independent, which has been shown to be true. As what was already said, 
the smaller the masking time is, the smaller the uncertainty is. However, instrumental constraints do 
not allow to rotate MicroSTAR too fast. For a given masking time, Fig. Ogives the optimal modulation 
period. For example, it shows that the set of parameters used until now (r = 10 min and Tm = 200 s) is 
"optimal" for the GLS approach, i.e. r = 10 min gives the minimum uncertainty for a masking time of 
200 s. 

6 Conclusion 

The Gravity Advanced Package, developed to improve orbit reconstruction of interplanetary probes in 
order to test General Relativity, relies on a technological progress with allows using an electrostatic 
accelerometer to make measurements with no bias. Indeed, the addition of a rotating platform allows 
modulating the non-gravitational acceleration while keeping the bias at low frequencies. The data acquired 
need to be processed in order to obtain the measurement with no bias. This data processing was the 
topic of this article. 

The first result obtained was conditions under which the bias is completely removed from the signal 
of interest. These conditions allowed designing calibration signals, i.e. a time-pattern for the rejection 
angle. Then the uncertainties on the unbiased non-gravitational acceleration were computed. It was 
shown that it is possible to recover the mean acceleration over each period of modulation and to have 
access to sinusoidal variations of this accelerations with some restriction on the pulsation of the signal. 
Finally, a method was presented to optimize the modulation period and the masking time so as to reach 
the minimum uncertainty. 

It has been shown that several parameters influence the precision on the post-processed quantities: the 
modulation time, the masking time and the integration time. It is possible to choose a set of parameters, 
which are technologically speaking reasonable, leading to precision below 1 pm s~ 2 on mean quantities as 
well as on the amplitude of sinusoidal variations. This precision is expected to improve orbit reconstruction 
significantly. 
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